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A theoretical analysis of signal fading using an energy density antenna 
is developed and compared with that from an isotropic antenna. The 
energy density antenna provides a signal proportional to the energy density 
of the mobile radio field. The number of crossings that the signal makes of 
a given signal level and the average duration of fades below a given signal 
level have been derived theoretically for these two cases using a simple sta- 
tistical model. Comparing the number of level crossings of the electric field 
with that of the energy density, it is shown that the energy density fades less 
frequently than the electric field by at least a factor of two. The average 
duration of fades of the electric field is greater than that of the energy density 
only for lower signal levels. These results are in reasonable agreement with 
experimental measurements. 

I. INTRODUCTION 

The study of signal fading appears to be very important to mobile 
radio systems. When a steady sine wave is sent out from a fixed station, 
the signal received by a mobile receiver in motion fluctuates, or, in 
radio jargon, fades. The received signal fluctuates more rapidly as both 
the frequency of the transmitted wave and the speed of the mobile 
radio increase. For a field received by a moving isotropic antenna, the 
maximum fading frequency f d , as Ossanna 1 has pointed out, is f d = 
2V/\, where V is the speed of the mobile radio and X is the wavelength 
of the steady sine wave. For instance, at 836 MHz and a speed of 15 
miles/hr, the signal fades at a rate of about 40 times every second and 
is a serious disturbance to the mobile radio communication. 

There have been many investigations of the fading problem. Aikens 
and Lacy 2 made a test using 450-MHz transmission to a mobile receiver 
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in an urban area. Bullington 3 investigated radio propagation variation 
at VHF and UHF. Young 4 pointed out that for the test samples of signal 
strength taken over a small area, the amplitude follows a Rayleigh 
distribution to a fair approximation. Also S. 0. Rice 5 pointed out that 
the fluctuations of a received radio signal have the same behavior as 
the envelope of a narrow-band Gaussian noise. Recently, Ossanna 1 
measured the power spectra of a mobile radio fading signal. They all 
treated the signal as obtained from an isotropic antenna. 

In this paper, a theoretical analysis of fading using an energy density 
antenna is developed and compared with that from an isotropic antenna. 
The concept of using the energy density antenna to reduce the effect 
of signal fading was suggested by J. R. Pierce. 7 It will be discussed in 
detail later. The number of crossings n(¥) that the signal makes of a 
given signal level ^, and the average duration of fades t(^) below a 
given signal level V, have been derived theoretically from a statistical 
model using Gaussian random amplitudes and equal angles of arrival 
of an infinite number of incoming waves. The two statistical factors, n 
and t, first expressed by Rice, 8 can describe the property of individual 
signal fading very well. In this paper, n and t for the isotropic antenna 
will be compared with the values for the energy density antenna. These 
theoretical results also will be compared with the experimental data. 

II. THE METHOD OF OBTAINING THE EXPECTED NUMBER OF LEVEL CROSS- 
INGS AND AVERAGE DURATION OF FADES 

From Kac's and Rice's 8 paper, a brief derivation of the expected 
number of level crossings n(*) of a given signal level ^ and average 
duration of fades below a given signal level ^ is as follows. We assume a 
random function \p which is statistically stationary in time, and for 
which the joint probability density function of \f/ and its slope \j/ is 
p(}J/, ip). Any given slope ip can be obtained by 

* = ijL , (i) 

T 

where t is the time required for a change of ordinate #, as shown in 
Fig. 1. The expected number of crossings of a random function \p in 
the interval (¥, SF — dip) for a given slope $> in time dt is 

the expected amount of time spent in the interval 
dip for a given \f/ in time dt 



the time required to cross once for a given \p in 
the interval d\p 



at >A-* 
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Fig. 1 — The notation used in obtaining the expected number of level cross- 
ings n(^) and the average duration of fades t(*). 



_ E(t) _ p(xl/, ^)d\pd^dt 
t d_\p 

4 



= M*. fidjdt. (2) 



u v--* 



The expected number of crossings for a given \f/ in time T is 
/ fo(V, >p)d\pdt = fat?, fydjsT. 

Jo 

The total expected number of upward crossings in time T is 

Jo 

The total expected number of crossings per second is 



(o; 



(4) 



(5) 



Since the expected number of crossings at a particular level Sl> per second 
can also be stated as 



n(¥) = 



the expected amount of time where the 

functio n \p is below level jg in one second 
the average duration of fades below level ^ 



hence, the average duration of fades below level ^ is 

PU < go 



t(*) = 



n(*) 



(6) 



(7) 
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Hence, the results will be derived from the joint probability density 
function p(\f/, r/-), and the problem is to derive this probability density 
function for the various signals. 

III. THE EXPECTED NUMBER OF LEVEL CROSSINGS AND THE AVERAGE 
DURATION OF FADES FOR A VERTICALLY POLARIZED WAVE 

In order to obtain the expected number of level crossings of a given 
signal level R and the average duration of fades below a given signal 
level R for the three field components of a vertically polarized wave, 
first we need to specify the forms of the three field components. Then a 
statistical model of the field components is assumed. From such a 
model, we find the joint probability density functions of amplitude R 
and its slope R for the three field components. Finally, we use (5) and 
(7) to obtain the result for each field component. 

Following Gilbert 10 a vertically polarized plane wave E, traveling 
in a direction u in the (x, y) plane is assumed. The three field compo- 
nents referenced to a receiver moving with velocity vector V can be 
written 

E z = e z = A u exp (-j/3u-v7) exp (jut) volt/m 

H x = v(h x amp/m) = A u s'm 6 U exp (-;'/3u-V7) exp (jut) volt/m 

H„ = v(K amp/m) = — A u cos 6 U exp (— jfii-Vt) exp (jut) volt/m, 

where /3 is a wave number and A u is a complex amplitude of an electric 
wave propagating at a direction u. u is a unit vector related to an angle 
U between the positive z-axis and the unit vector itself. 77 is free-space 
wave impedance. The time variation exp jut can be dropped out of 
three field components for simplifying the derivation. Moreover, from 
now on, we will treat the units of all three components E, , H x , and 
H y in volt/?n which will also simplify the calculation. 

When N vertical polarized waves coming from AT directions are re- 
ceived by an isotropic antenna of the mobile radio, the three components 
become 

N N 

E.= 22 A « ex P (-#u-V0 = J2 A u exp [-jpVt cos (6 U - a)] (8) 

u =l u=l 

N 

H x = 22 A u sin 6 U exp (-j/3u- Vt) 

u-l 

If 
= 22 A u sin 6 U exp [-jpVt cos (0 U - a)] (9) 
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H y = £ - A u cos U exp (-J/3U-V/) 

= - 2 4„ cos 0„ exp [-jfiVt cos (0„ - a)], (10) 

where U is the angle between the positive x-axis and the direction of 
wth wave u, and ^ U ^ 2ir. a is the angle between the rc-axis and the 
velocity V, and ^ a ^ 2t. Both U and a are shown in Fig. 2. 

In this paper, a statistical model is used as follows: The complex 
amplitude A u can be separated into a real and an imaginary part A u = 
#u + i'Su , hence N incoming waves have N real values of R u and S u . 
We suppose all those 2N real values are Gaussian independent variables 
with mean zero and variance one. Also, we assume the N waves have 
uniform angular distribution, i.e., the fcth wave u k has an angle of arrival 
U = 2irk/N. Moreover, in this paper an infinite number of multiply 
reflected waves (N — * «• ) are assumed for finding the expected number 
of level crossings n(R) of a given signal level R, and the average duration 
t(R) of fades below a given signal amplitude R. 

3.1 Finding the Values of n(R) and t(R) from the E, Field 

First of all, we need to obtain the joint probability density function 
of signal amplitude R and its slope R for the electric field component 
E t using the statistical model we mentioned previously. We start from 
(8). The alternate form of (S) can be written as 

N 

E. = Z (Ru + A) [cos {fiVt cos (0 U - a)} - j sin [fiVt cos (0« - a)}]. 

(ID 




V = VEHICLE VELOCITY 



U= DIRECTION OF PROPAGATION OF 
A RANDOM FIELD COMPONENT 



Fig. 2 — The coordinate system. 
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Equation (11) can be separated into real and imaginary parts 

E. = X, + jY> . (12) 

The real part of E z is 

Xx - T, (B- cos <p u + S u sin <p u ) (13) 

and the imaginary part of E, is 

Y l = £ (& cos *> u - fl„ sin <p u ) , (14) 

u = l 

where 

^ u = pVt cos(0 u - a). (15) 

We assume that all N values of R and S in (13) and (14) are time in- 
dependent. Then the derivatives with respect to time of (13) and (14) 
are 

AT 

-X\ = pv 23 (— R* sin ft + s « cos ft) cos (*» ~~ a ) v16 ) 



II 1 



F, = /37 2 C" 5 " sin ft ~~ ^ u cos ft) cos ^ u ~ ")■ ^ 17 ) 

II 1 

The mean values, variances, and covariances of X, , F, , X x , and F, are 
mi = (X t ) = (Y l ) = (X,) = (Y,) = 
Mll = (XI) = <F 2 > = N for any tf (18) 

rf, = (A1> - («> = (W f for tf ^ 3 (19) 

and 

(Z.F,) - (XX) = (Y^) = (Y l f l ) = (X.F,) = (X,F0 = 0. 

The above results are shown in the Appendix. 

From the central limit theorem, it follows that Xi , Y\ , X i , and Y 1 
are four independent random variables which are distributed normally 
as the value N approaches infinity. The probability density function 
of four independent real random variables X x , Y x , X, , and F, is 11 

p(X, ,Y ll X 1 , F,) 

1 / 1 ( X] + F; , Al + II 

= (2^i^ exp r 2 r^r - + ~^tt 
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where the determinant |m| of the covariance matrix is 

We may introduce the concept of the envelope 
E, = X l + jY 1 =rjju. 

The quantity r, is the envelope and 77 „ is the phase, both of which are 
slowly varying functions of the time. Then, 

X, = r. cos 77. ; Y x = r e sin 77, (21a) 

Xi = r. cos 77. — r.T). sin 77, ; F, = r e sin 77, + r.17, cos 77. . (21b) 

The Jacobian of the transformation from (X! , Y^ , X t , Y\)-space to 
(r. , v. , r. , 77.)-space is 12 \ J \ = r*. 

Therefore, the change of variables gives the probability density the 
form 

p(X 1 , Y 1 , X, , F.) = r*q(r t ,v.,K, v.) 

= p(r. ,V,,K, v.) 

_d__ exp (_I(i + %Lil)|, (22) 



(2?r) I /* I I 2 Vn M11 

where g(r,, , -q, , r, , 17,) is the density obtained on substituting for 
Xi , Fi , etc., their values in r. , 77. , etc., obtained from (21a) and its 
time derivative (21b). To obtain p{r. , f«), the probability density of 
the envelope and its rate of change, we must integrate over 77. and 77, , 
the range of which are, respectively, (0 to 2r) and (—00 to + °° )• From 
(22) we obtain 

p( r . , f.) = -7==- exp {l (^ + 4)}- (23) 

It is observed that the expression on the right of (23) is independent 
of t. Hence, the expected number of level crossings n{R t ) at a given signal 
amplitude (r, = R.) can be obtained from (5) by using p{R. , r.) in (23). 



<R.) - ftpiR. ,t)dt = >/f d==exp {-&) 

Jo ^Mn V 2*7*11 X -M11/ 

Now the variance of r, is 

<r?>- (X 2 ) + <F 2 ) = 2 Mu = 2AT. 



(24) 
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R. = R./Vtf) = BJr. trmti = R./V2N. (25) 

Substituting the values of variances mh and n' u from (18) and (19) into 
(24), also applying the relations in (25), we obtain 



nfc) = -££=#. exp(-jBJ). 



(26) 



Equation (26) is plotted in Fig. 3 where the abscissa is R. in dB (20 log 
R e ) and the ordinate is ( VW|8V>(#.)- 

The average duration of fades t(R.) of E, can be obtained as follows: 
The probability that the envelope r e is less than a given amplitude level 
R t is 

P(r.(0 <R.) = [ ' V{r£t))dr.(t), (27 ', 

Jo 
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Fig. 3 — Comparison of the level crossing rate of the electric field with that of 
the energy density. 
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where p(r,(t)) is the probability density obtained from p(X x , Fi). By 
changing the variables, we obtain p(r e (t), r) e (t)) from p(X t , Y t ). Then 
we integrate p(r t (l), 77. (J) ) over from to 2x to obtain 13 



P(r t (0) = 7^ exp 



-f 1 • 



where /uu is obtained from (18). Substituting (28) into (27) wo get 

P(r e (l) < R.) = 1 - exp (-R;/2fi u ) = 1 - exp (-/?;). (29) 

The expected number of times per second that r e (t) passes upward (or 
downward) across the level R e is obtained from (24). The average dura- 
tion of fades during which r e (t) < R e may then be obtained by sub- 
stituting (24) and (29) into (7) 

= PMl^EA m E<m<lA . ^ i [exp m _ „ (30) 

n(it„) n(R e ) 0V R e 

which is shown in Fig. 4. 

3.2 Finding the Values of n(R hx ) and t(R hI ) from the H x Field 

Following the same steps as above, we are going to find the joint 
probability density function p(r hz , r hx ) of the envelope r hx and its slope 
r hx of the II x field component first. From (9) we find the real and imagi- 
nary parts of H x which are expressed in the Appendix. The means, 
variances, and covariances of four real Gaussian random variables 
shown in the Appendix are 

m 2 = (X 2 ) = (Y t ) = <X>> = <f 2 > = 
M22 = X\ = Y\ = — for any N 

M 2 ' 2 - Al= Y\ = | 03F) a 

• [cos 2 a + 3 sin 2 a] for N = 3 and N ^ 5 

(X 2 Y 2 ) = (X 2 X 2 ) = (X 2 Y 2 ) = (Y 2 X 2 ) = (Y 2 Y 2 ) = (X 2 Y 2 ) = 0. 

The probability density of the envelope of H x field and its rate of change 
pi^hx , r hx ) is then obtained by following the same procedure used in 
deriving p(r t , r,). 

p(r hI , h,) = —?^=r ~ ex P \ - o yr + ;W f • ( 31 ) 

V27T/Xo 2 M22 I. ^ V*22 M22/J 
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• t(* t ) THE AVERAGE DURATION 
OF FADES BELOW A GIVEN LEVEL 
$ t OF THE ENERGY DENSITY 

▲ t(Re) THE AVERAGE DURATION 
OF FADES BELOW A GIVEN LEVEL 
Re OF THE ELECTRIC FIELD 




-20 -15 -10 -5 5 10 

SIGNAL LEVEL IN dB WITH RESPECT TO RMS VALUE 

Fig. 4 — Comparison of the duration of fades of the electric field with that of 
the energy density. 

Hence, the expected number of level crossings n(R hx ) at a given signal 
amplitude (r hz = R hx ) can be obtained from (5) by using p(R hx , r hx ) 
in (31) 

«QU = jffcrffc , W*. - ^g (^=) oxp (-£)• (32) 

The variance of r hz is 

(rl) = (XI) + {Yl) = 2 M22 = N. 
Let 

Rhx Rhx Rhx 



Rhx = 



V(rl) 



r h x(rm.) 



Vn 



(33) 
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Substituting the values of the variances M22 and M22 into (32) and replac- 
ing AT by (r 2 hz ) we get 

BV /cos 2 a + 3 sin 2 a 



n(R hx ) = 



/2k 

BV 
V2x 



2 



' R hx exp —R\ 



-\/l — \ cos 2a R hx exp —Rl x . 



(34) 



Equation (34) is the same form as (26) except for a multiplying factor 
which is a function of a shown in Fig. 5. Hence, n(R hx ) is also a function 
of angle a. Thus, when the mobile is moving along the z-axis a = or 
ir, and 

BV 1 « 
V%r V2 



n(R hx ) = —7= —7= R hx exp -R 2 hl 



(34a) 



which is the minimum value of n(R hx ). When the mobile is moving on 
±y-axis a = ±ir/2, and 



n(R hx ) = 



BV 



- R hx exp — R 2 hx 



(34b) 



which is the maximum value of n(R hx ). The ratio of level crossings for 
these two cases is 

n(R hr )(<* = 0°, 180°) = 1 
n(R hx )(a = ±90°) VjT 



(34c) 
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Fig. 5 — The effect of the angle a on the ratio of level crossing rates of the 
electric field to the z-component of the magnetic field. 
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For a = ±45° and ±135°, 

n(R hx ) = -$j= £ te exp -Rl (34d) 

which is the same expression as (26) for R. . However, although (34d) 
and (26) are the same form, the magnitudes R e and R hx are different, 
since (r 2 e ) in (25) is equal to two times (r 2 hz ) in (33). 

The average duration of fades t(R hx ) of H x can be obtained without 
difficulty. It is easy to prove that the expression for the average length 
of the intervals during which r hx < R hx will have the same form as 
t(R.) in (30), except for a multiplying factor that depends on a, as 
follows: 

tCSj = ^ -t; \ , J" [exp (Rl) - 1] . (35) 

pv V 1 — i cos 2a Khz 

When a = 0° and 180°: 

t(R hx ) = -^ V2~ [exp (Rl) - 1] (35a) 



which is the maximum value of t(R hz ), and when a = ±90°: 



(35b) 



which is the minimum value of t(R hx ). The ratio of level crossings for 
these two cases is 

t(g ta )(«-0°,180°) = ^ (35c) 

t(5 to )(« = ±90°) 

which is the inverse of (34c). This tells us that when n(R hx ) reaches a 
maximum value, the average duration of fades reaches a minimum value 
and vice versa. 

For a = ±45°, ±135° 

t(R hz ) -%£■ [exp (Rl) - 1] (35d) 

pv n hx 

which is the same form as the expression for t($„) in (30). We may say 
at these angles a = ±45° and ±135°, the E field and the H x field have 
the same average duration of fades below the level Ti, = R hx . 
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3.3 Finding the Values of n(R hu ) and t(R hy ) from the H y Field 

Similarly, we obtain n{R hy ), the expected number of level crossings of 
a given signal level R hy from the H y field. It is very easy to see that 
n(R hz ) and n(R hu ) are the same forms of distribution as expressed in 
(34), except for the multiplying factor that depends on a. The average 
duration of fades t(R hy ) during which r hy (t) < R hy is also of the same 
form as t(R hx ) in (35) except for the multiplying factor depending on a. 
The variances /n 33 and n' 33 are given in the Appendix. We may thus write 
directly 

/3T 



n(R hv ) = - £ r= Vl + \ cos 2a R hy exp -Rl„ (36) 



where 



1 /? irr /, ■ \ o h [exp ^u " 1L (37) 

py Vl + h cos 2a R hy 



R* 



It is obvious that 



when R hy = R hx . 



(fly) rhvfrm.) 



(Rhvja-o'.iBo' — n{R>,x)a-±00' 

n(Ri, y )a = ±90° = ft\Rhi)a-0°.lSO'' 



IV. THE EXPECTED NUMBER OF LEVEL CROSSINGS AND THE AVERAGE 
DURATION OF FADES OF THE SIGNAL FROM AN ENERGY DENSITY 

ANTENNA 

J. R. Pierce 7 has suggested utilization of the energy density concept 
as a possible means for reducing the signal fading in mobile radio. If 
we pick up the electric field e and the magnetic field h in free space and 
amplify the two fields by their appropriate relative gains, square and 
add these two fields, we obtain a signal proportional to electromagnetic 
energy density 

W = \{ee 2 + y.h 2 ), (38) 

where e is dielectric constant, and m is permeability. This idea can be 
realized by using a special antenna which receives three field compo- 
nents e t , h z , and h y simultaneously. The three signals enter separate 
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square-law detectors, and the three detector outputs are added to obtain 
the energy density, 

w - }(■ M + /* W* + m W). 

We may express IF in a different form 



2 



- - e p.| 8 (voltVm 8 ) + \H X \ 2 (voltVw 8 ) + [ffj (volt 2 /™ 2 )] 



= s [^.(volt'/w 2 )] - *£ Joules/m 3 . (39) 



2 
We define \f/, as a normalized energy density 

*, = |£,| 2 + |ff,| 2 + |tf,,| 2 volt 2 /m 2 

= fc + ft, + *» W 

= (X 2 + F 2 ) + (X 2 + Y\) + (X 2 + F 2 ). (41) 

Gilbert 10 has done some work on finding power spectra in energy re- 
ception for mobile radio. His work provides very useful background 
for this paper. 

In this section, we are attempting to derive the number of crossings 
n(¥.) at a given level of signal magnitude V, using (5) in Section II. 
First of all, we need to find the joint probability density function 
Pitt Ad of signal f t and its slope ^, . Since f, is a function of 
(X, , F, , X 2 , Y 2 , X 3 , F 8 ), and & is assumed to be a function of 
tt. , i'hx , thy), we will nn d out that the variables (X t , Fj , X 2 , F 2 , 
X a , F :) ) and 0/-. , \p hz , W are two independent Gaussian variable 
groups. Then, 

p[^ t (X 1 ,Y l ,X 2 ,Y 2 ,X 3 , F 8 ), &(& , K , K)] 

= p[^(X, , Y l , X a , Y a , X, , Fa)] X vlM4>. . fc- . Wl 

- p(h)vM- <&> 

A brief sketch of the method of finding p(tf-,) and p(^«) is discussed below. 

4.1 To get pO/0 

Since we know from (41) that 

*, = (X 2 + Yl) + (X 2 + F 2 ) + (XI + Yl) 
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and since X, , F, , • • • are independent Gaussian variables, it is easy 
to get p(\p,) from the Fourier transform of the characteristic function 
M+ t (v), where M+ t (v) = M x M Y M x M Ya , ••• etc., and we can get 
these il/'s very easily. 

4.2 To gel p(\p t ) 

All the terms in the summations of equations (79), (80), and (81) 
which represent \j/ e , \p hx and \p hv , respectively, in the Appendix are 
statistically independent. Then by the central limit theorem these 
three variables \p e , ij/ hx , and \p hu are Gaussian distributed. Hence, the 
joint probability density function p(\js e , \J/ hx , ^ hv ) can be established. 
Since \j/ t = $, + $ hI + ^ hy we can get p(\p,) from the Fourier transform 
of the characteristic function 71/ ; ( , but il/^ ( must now be obtained from 
the general definition 



M h (v) = E[e"*>] - JfJ V'Stf. , K , KWMM h 



(43) 



since we have no simple way of getting Mf t , M^ km and M^ ht separately. 
Let us introduce a new variable e which can be any one of the above 
Gaussian random variables. It has a zero mean and variance /x. Then 
the probability density function of the square e is 1 ' 1 

for 7 > 0. The characteristic function corresponding to this probability 
density is 



M y (jv) = C e ,ty p(y)dy = (1 - j2^ . 
Jo 



(45) 



From the Appendix we know all six variables X x , Y x , X, , Y 2 , X 3 and 
F 3 are independent Gaussian variables. It is not hard to see that the 
X] , Y\ , XI , Y\ , X\ , Y\ are independent variables by obtaining 
p{X\ , Y\ , X\ , Y\ , X\ , YD from the Jacobian of the transformation 12 
of p(Xi , F, , X 2 , F 2 , A r 3 , F 3 ). Then X\ and F 2 have the same char- 
acteristic function (1 - 2iMn*>)~*. Also X\ , Y\ , X 2 , and F 2 have the 
same characteristic function (1 - 2;n 22 v)~ i . Thus, by the addition theo- 
rem 15 the sum \p t which is defined in (41) has the characteristic function 

**** - (i - 2m J { i - 2 iM ,f m 
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Then the probability density function p(^«) can be obtained from the 
Fourier transformation of the characteristic function 

a, r ( m / h w , *i exp (a 

= 2( M . - m„) 2 exp rw " cxp rw. + w*. - *o 

(47) 

The joint probability density p(& , ^ , &») has been derived in 

the Appendix 

p(^« . ^*x . thy) 

= 7^-Tn-T~TI CX P ( _ KI A " I ^ + I A 22 I & + I A 83 | &„ 
(27r)* I A | a 

+ 2 | A 12 | ^ Aj! + 2 | A 13 | t.K + 2 | A 23 | ^Wl, (48) 

where [A] is the co variance matrix of $, , \{/ hz , j/- Av and the |A» m | are the 
cof actors of |A|, given in the Appendix. From (40), ^, is the sum of 
the three random variables 0. , iA Aj , and ^», . Then the characteristic 
function for \j/ t is 

CO 

= f/j exp {jv(\p. + lA*. + W)p(& , \ta . ^hv)d^ e d^hzdyl/ h y . 

(49) 

The details of this computation are given in the Appendix with the 
result (92) 

MiJjo) = exp {-\ P [v 2 \, (50) 

where 

P ', = Pl'l + P22 + p^ 3 + 2 P [ 2 + 2 P ' 13 . (51) 

The probability density of total yj/ t is then 

*« ■ s £ *'.«•-*'* - vb exp (~i '•)• (52) 

The joint probability density p(^, , ^) of iA, and ^, can now be ob- 
tained by substituting (47) and (52) into (42). 
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The expected number of crossings n(^ t ) at a given signal level \f/ t = 
tyt made by the total energy density signal in one second can be ob- 
tained from (5) 

n(*0 = f t,p(% , tM t = p(%) f i,p(h)dh 

Jo Jo 

•V, exp -^,/2m 2 2 



Mm 



27T L2(M22 - Ml.)' 



exp 



+ 



4/J 22 (m22 — Mil) 



(53) 



where p' t is given by (51) and Mn = Af, M22 = N/2 as shown in the Ap- 
pendix. 
Also we know from Gilbert 10 

<tf) = 22N 2 . 

In addition, p{ = p' u + p 22 + p^ 3 + 2p( 2 4- 2rf 3 = 2N 2 (0V)' 2 and 
letting 



*, = 



*! 



*, 



we can simplify (53) as follows: 
n($ t ) = -^=J2V2 exp (-^^^j - exp (- V22$,) 



- 2Vli$ t exp (- \ / 22% l )>- 



(54) 



Equation (54) is a distribution which is independent of the angle 
a. When $r t = 1, it means that ~$r t is equal to its rms value. Equation 
(54) then becomes 



«(#)*., = 



Also let/?, = 1 in (26): 



ra(#.)5.=i = 



0F 



07 



X 0.1839. 



X 0.3678. 



It is shown that the expected number of crossings of the total energy 
density is one half the expected number of crossings of the envelope of 
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E field at their rras levels ($, = R 2 a ). In other words, the energy density 
fades less frequently than the E field by at least a factor of 2. 

for 4f, = #* with respect to their rms values. The theoretical values 
of n (^) and n(R.) are in Fig. 3. 

The average duration of fades below a given level ^r t is given by (7) 

«w - ^^ ■ « 

where P(yp,(t) < ^,) is the probability function obtained from p(\p t (t)). 
P(+ t (t) <*<)=[' p(Mf))dM) 

Jo 

V22 



= 1-4 exp (— ^ *.) + (3 + V22*0 exp (- V22*,) , (56) 
where 

*.- *' 



Substituting (54) and (56) into (55), we obtain the average duration 
of a fade below a given level ty t : 

1 - 2 exp (-"^ $t +) exp (- V22^«) 



-1 



2[exp f — ^y^ *J - exp (- V22&«)]- V22* ( exp (- V22*<) 

(57) 
When ^ = 1, (57) becomes 

t(*0 = ^ x 3 - 74 - 

When #, = 1, (30) becomes 

t(R e ) = ^r X 1.7183. 

It is shown that the average duration of fades of the energy density 
below a level ty, is larger than the average duration of fades of the E 
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field below a level R. where % e = R] = 1. The curves of t(#,) and 
t(R e ) have been plotted in Fig. 4 for comparison. 

V. DISCUSSION OF THE THEORETICAL RESULTS 

From the above derivation we know that \f/ t and \p t are two independent 
variables as shown by (42): 

v(ti . h) = pMpM, 

and therefore (5) can be written as follows: 

Jo 

= p(*«){lM (58) 

where {i/-,} represents an integral. Equation (58) simply shows that 
the expected number of crossings n(^ t ) at given level ^ can be obtained 
from the probability density of level ^, times the integral {^,} . The 
average duration of fades, then, turns out to be 

W n(¥ ( ) 

We emphasize that (58) and (59) are valid only when \f/, and \p t are two 
independent variables. 

The two curves, n(R.) and n($ t ), are plotted in Fig. 3, normalized 
by the common factor VW/SF. Both curves are plotted as functions of 
the signal level normalized to their own rms values. The value of n(R„) 
is, as shown, always higher than the value of n($ t ) for any signal level. 
From these two curves, it may be said that the fading of the energy 
density is less frequent than the fading of the envelope of the electric 
field. The maximum expected numbers of crossings of both n(R„) and 
n(&t) are at the -3 dB level, which means for signal level at \/\/2 and 
i of their rms values, respectively, we will count the most fades. The 
curve of n($ t ) has dropped faster on both sides of dB than the curve 
of niR,), which means that the range of the signal amplitude \p t is 
less than the range of the signal amplitude r e . 

The average duration of the signal below a given amplitude level 
is another way of looking at the fading problem. Fig. 4 shows that the 
average duration of fades of the energy density t(^,) is always larger 
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than the average duration of fades of the electric field t(R t ) when the 
given signal levels are above — 3 dB with respect to their rms values 
(*. = R] > -3 dB). The value of t(*,) is less than t(R.) when the 
given signal levels are more than 3 dB below the rms values ^ t = R 2 e < 
-3 dB. When the given signal level is at -3 dB (*, - j£J - -3 dB), 
the average duration of fades t( — 3 dB) of both the energy density and 
the electric field are the same. 

VI. COMPARISON OF THE THEORETICAL PREDICTION WITH THE EXPERI- 
MENTS 

The three field components E, H x , and H v have been received by a 
special antenna 0,7 mounted on a mobile van moving at a speed of 15 
mile/hr. All the figures shown in this section were taken on Common- 
wealth Avenue, New Providence, New Jersey, from a transmitting 



ta 



I* 




^-15 -10 -5 5 10 

Re IN dB WITH RESPECT TO RMS VALUE 



Fig. 6 — Comparison of the predicted level crossing rates to the observed rates 
for the electric field on Commonwealth Avenue, New Providence, New Jersey. 
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-15 -10 -5 5 10 

Re IN dB WITH RESPECT TO RMS VALUE 



Fig. 7 — Comparison of the predicted average duration of fades to the observed 
average duration of fades for the electric field on Commonwealth Avenue, New- 
Pro vidence, New Jersey. 



antenna at 836 MHz at Bell Laboratories, Murray Hill. After adjusting 
the appropriate relative gains of the three fields, the energy density 
can be obtained by squaring and summing these three fields by computer 



* ( = | E | 2 + | H z | 2 4- | H v | : 



volt 2 /™ 2 



Since the distance between the transmitting antenna and the mobile 
unit is relatively short, the angle swept out by the radius vector from 
the base station to the mobile unit varies considerably over a typical 
length of run. To reduce the variation of this angle the data for the 
entire run were cut into sections 8 seconds long, corresponding to 175 
feet of travel, for computer processing. Each section, either the en- 
velope r. of the E field or the energy density \J/ t , was used to obtain the 
number of level crossings n and the average duration of fades t by com- 
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-20 ^-15 -10 -5 5 10 

+ t IN dB WITH RESPECT TO RMS VALUE 

Fig. 8 — Comparison of the predicted level crossing rates to the observed rates 
r the energy density on Commonwealth Avenue, New Providence, New Jersey. 



puter program. However, since the experimental curves of n and t were 
almost all alike for all sections, we used only one for comparison with 
the theoretical curve. 

Fig. 6 shows a comparison of the curves of the expected number of 
crossings n(R.) at any level R a for both experiment and theory. The 
shape of the experimental curve is in fairly good agreement with the 
theoretical curve. Since the receiving antenna on Commonwealth 
Avenue is in line of sight with the transmitting antenna at Bell Lab- 
oratories, a small direct wave component may be introduced. This 
small direct wave component is not considered in our theoretical analy- 
sis, hence the values n(R e ) from the experiments should be less than the 
theoretical results as we would predict. 

Fig. 7 shows a comparison of the curves of the average duration of 
fades t(R e ) for both experiment and theory. They are quite alike. Since 
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a small direct wave component does exist, the average duration of fades 
for the experimental data should be higher than the theoretical results. 

Fig. 8 shows a comparison of the curves of the expected number of 
crossings n(^«) at any level ty t for both experiment and theory. The 
shape of the experimental curve is veiy much like the theoretical curve. 
It shows that the theoretical model used in this paper is quite acceptable. 

Fig. 9 shows a comparison of the curves of the average duration of 
fades t(#,) for both experiment and theoiy. The difference between the 
experimental curve and the theoretical curve may be caused by the 
small direct wave component. A small direct wave component intro- 
duced into our theoretical model may cause a little higher average dura- 
tion of fades than it might expect, but does not affect the number of 
level crossings. 
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Fig. 9 — Comparison of the predicted average duration of fades to the observed 
average duration of fades for the energy density on Commonwealth Avenue, New- 
Providence, New Jersey. 
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VII. CONCLUSIONS 

Comparing the expected number of level crossings and the average 
durations of fades of the energy density with that of the E field, we see 
that the fading of the energy density is much less severe than the fading 
of the envelope of E field. 

Referring to Fig. 5, which shows the fading rate related to the orienta- 
tion of the energy density antenna and the direction of vehicle motion, 
we see that when the two orthogonal loops are at 45° to the direction of 
motion, the fades of all three field components are the same. When one 
loop is lined up with the direction of motion and the other normal to 
it, the H field component received from the loop normal to the motion 
has less fading than either of the other two field components. 

The expected number of crossings/second of fades at a given signal 
level, n, for both R. and •&, is proportional to the carrier frequency / e 
and the mobile speed V, as shown in (26) and (54). They have the com- 
mon factor, (/37/ V2» = V2^(Vf c /c), where c is the velocity of light. 
Hence, if either V or / e goes higher, n becomes greater. 

The average duration of fades, t, is inversely proportional to the 
carrier frequency / c and mobile speed V, as shown in (30) and (57). 
Hence, if either V or /„ goes higher, t becomes smaller. 

The foregoing theoretical analysis is based on a Gaussian model and 
does not include a direct wave component. Even so, this analysis is 
compared with the experimental results in Section VI with fairly good 
agreement. 
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APPENDIX 

A.l Finding the mean values, variances and covariances oj nine variables 

From (8), (9), and (10) we may express in the following forms 
E z = X x + jYt 
H x = X 2 + jY, 
H v = X 3 + jY a , 



where 



also 
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N 

X, = £ R« cos *» + s » sin *« (60) 

u=l 

N 

F, = X & cos <*> u - R u sin $ u (61) 

AT 

X 2 = J^ (^u cos 3> u + S u sin $„) sin 6 U (62) 

u = l 
Y 

F 2 = 2 (S„ cos $ u - # u sin *„) sin d u (63) 

u = l 

AT 

X 3 = - X) (#» cos * u + 5 U sin *„) cos 6 U (64) 

N 

F, = - £ {S u cos 3> u - i2„ sin $„) cos ft, (65) 

u=l 

* u = (87/ cos (ft, - a) (66) 



and the angles ft, and a are shown in Fig. 2. The time derivatives of 
X t , Y r , X 2 , Y 2 , X 3 , and F 3 are 

^i = fiV £ (-72. sin $ u + £ u cos *J cos (d u - a) (67) 

7, = 07 X (-& sin * u - fl„ cos *„) cos (ft, - a) (68) 

& - 07 2 ( _ #« sin *■ + & cos *«) sin *« cos (ft, - a) (69) 
F 2 - /37 £ (-& sin <J>„ - #„ cos *«) sin ft, cos (ft, - a) (70) 

X 3 = -07 2 (--R u sin * u + £„ cos $„) cos ft, cos (ft, - a) (71) 

F 3 = -07 2 (-& sin <!>„ - /?„ cos *„) cos ft, cos (ft, - a). (72) 

The mean values of all above random variables are zero (X t ) = (F,) = 
(X 2 ) = (F 2 > = (X 3 ) = (F 3 > = (Z,) = (t t ) = (X 2 ) = <F 2 > = (X 3 ) = 
(F 3 ) = 0. The variances of all above random variables are 
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Mn - (XI) = (Y 2 i ) = N 
M22 = (Xl) = (Yl) = -r 

n 33 = (X 2 3 ) = (Yl) = — 



rf, = <X 2 > = <y 2 > = |(/37y 



(73) 
for any JV (74) 

(75) 
foriV^3 (76) 



M*, = <X 2 2 ) = <F 2 2 ) = ? (/37) 2 [cos 2 a + 3 sin 2 a] 



Ma 



33 - <i!> = (Yl) = ? (/3F) 2 [3 cos 2 a + sin 2 a] 



- for , 



iV = 3 (77) 
liV ^ 5.(78) 



Remark: The values of Mn for iV ^ 3 is derived as follows: The sum- 
mations of sine and cosine functions can be expressed 



22 sin kx = 



cos - — cos (2JV + 1) o 

_ . x 
2 sin g 



yi cos kx = 



sin (2iV + 1) | - sin § 



2 sin | 



Then in (67) 



2ttm 
iV 



and 



2 cos 2 (0„ - a) = g Z [l + cos (u -^ - 2a] 



iV , cos 2a 



= 2 + 



47r , sin 2a 



Ar 



Z17T . sui ^a v> • tvr 
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But, 



,, sin ^ + |)- sin | 

£ C0Stt AT = ~7Y, = ° 



2 sin 



N 



and 



A . 4. COS ^~ COS l 47r + ^ 

2>n«3jr- — -^ 

2sin iv 

Thus, the average value of X, or 7? is 



•for N ^ 3. 



- 



(f> s °(- 



2 /2ttu 
2^ COS [-JT- - a 



= % for N £ 3. 



Following the same derivation, we obtain the valid range of AT for n' 22 

and ^33 . Later on we will obtain the range of N for p' n , p 22 , p 33 . P12 , P13 

and P23 in the same way. 

QED. 

Now we are going to find the relations between all the six variables 
X x , X 2 , X 3 , F, , F 2 , F 3 and their time derivatives. Since we know 
if two variables a and b are Gaussian, and also uncorrelated, (ab) = 0, 
then a and b are independent. 16 Therefore, the covariance of X x , F, , Xi , 
and Y x are 

<X,F,> = (XX) = (X 1 f l ) = (Y 1 X 1 ) = (Y 1 Y 1 ) = (X.F.) = 

hence, the four variables X, , Yi , Jti , and F, are statistically independ- 
ent. The covariances of X 2 , F 2 , X 2 , and F, are 

(X 2 Y 2 ) = (X 2 X 2 ) = (X 2 Y 2 ) - <F 2 X 2 > = (Y 2 Y 2 ) - <X 2 F 2 > = 

hence, the four variables X 2 ,Y 2 ,X 2 , and F 2 are statistically independ- 
ent. The covariances of X 3 , F 3 , X 3 , and F 3 are 

<X 3 F 3 > = (X 3 X 3 ) = (X 3 Y 3 ) = <F 3 X 3 > = <F 3 F 3 > = <X 3 F 3 ) = 

hence, the four variables X 3 , Y 3 , X 3 , and F 3 are statistically independ- 
ent. 

Also we may show that 

(X m F„) = for all m and n 

(X m X n ) = (Y m Y„) = for m 5* n, 
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where 



in 



n 



= 1,2,3, 



hence the six variables X x , X 2 , X 3 , F, , Y 2 , Y 3 are independent. The 
rate of change of energy densities of three fields E, , H z , and H y are 

= ^E [-(R U R V + S U S V ) sin (#„ - *.) 
+ (S U R, - RuS v ) cos(<S> u - *,)] X [cos(0 u - a) - cos(0, - «)], (79) 

= V Z [-CBJB. + &&) sin (<*>„ - $.) 

+ (S B fl, - R U S W ) cos (*„ - *.)] 
•sin 6 U sin 0,[cos (0„ - a) - cos (0„ —a)], (80) 

= PV 2 [-(«.«. + &<S.) sin (*. - *.) 

+ (S U R, - RuS v ) cos (<*>„ - $,)] 

•cos U cos r [cos (0 U — a) — cos (0, — a)]. (81) 

The only terms that exist in (79), (80), (81) are those for which u 5* v. 
There are N(N — l)/2 different terms which are all statistically in- 
dependent in (79), (80), and (81). Hence, by the central limit theorem, 
^. , ^h U > an d \p hx are Gaussian random variables. 
The variance of \j/, , \j/ hx and \j/ hv are 



Pi'i = <#> = (fiV) 2 *N(N - 1) 

p 2 ' 2 = <&> = (W ^JLZ 1) [cos * a + 3 sin 2 a] 

PL - <&) = (0V) 2 N(N ~ l) [3 cos 2 « + sin 2 a] 



for N ^3 (82) 



- for 



\N = 3 
AT > 5. 



(83) 



(84) 
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The covarianocs of \j/ e , yj/ hI , and \p hv are 
Pi'i = Pa. - (Mkz) = -WV)' 2 2N(N - l)shr a 
pU = pL = <M».> = ~(0V)*2mN - cos 2 a 
p 23 = pi. - <**.&,) ^ 0. 



for JV ^ 3 
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(85) 
(80) 
(87) 



It is very easy to show from (60) to (66) and (79) to (81) that the covari- 
ances of the variables between two groups (\p„ , fa, , \p hv ) and (X, , F, , 
X 2 , Y 2 , X 3 , F 8 ) are zero. We may write 



d m X n ) = {i m Y n ) = for all n and 



m 



(88) 



m 



= 1, 2,3 



hence (^. , ip hx , fay) and (X, , Fi , X 2 , Y 2 , X 3 , Y 3 ) are two independent 
variable groups. 1 

A.2 Derivation of Mj,,(J v ) * n (50) 

The mean values of all three Gaussian random variables i/-, , \j/ hz , 
and yj/ hy we observed from (79) to (81) are zero. Also (87) gives the 
covariance {^hx^hv) = 0. The joint probability density function of three 
variables ^„ , \l/ hx , ^ hv can be obtained 11 

p(& , K , K) = ( 2t )* I A i* exp I 2~[a~| ( ' Au ' ^ + ' A " ' ^*" 
+ I A M | & + 2 | A 12 | M ta + 2 | A 13 | *.&, 

+ 2 | A a3 | ^„„) j , 

where 

[A] is a covariance matrix, and 

|A m „| is a cofactor of p mn in the covariance matrix [A] 



[A] - 



Pll P\2 Pl3 

p[ 2 pi. o 

LPi3 paaJ 



(89) 
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|A| = determinant of [A] = p' 33 (p'nP22 ~ Pit) ~ PnPu 
A„ 

A 9 , 



A 33 
A 12 
A 13 
A 23 



— P22P33 



— PllP33 ~~ Pl3 
= PllP22 Pl2 



= I A 2 l I = — P12P33 
== I A 3] I = P22P13 

= I A 32 I = Pi' 3 Pi2 • 

From (38), r/- t is the sum of three Gaussian random variables \j/ e , \^x , ^a„ .' 

4*1 = Ipt + ^hx + thy • 

The characteristic function for \J/ t is then 
M h (jv) = £[exp {jv(+ e + K + &,}] 

/ / exp {jv(\jy e + $ hx + 4s hv )\v(^t - '/'ax , ^hy)diMh Z d^ hu 

•00 **— 00 *'— 00 

(2tt) j I A |* J_ M J_ M 
•exp {-§"]V[ ( ' Al1 I ^ + I A " I 4L + 2 I A» I &&,)} 

x £^M-^(*+ * |A "'fc ' i i,, '* J *-)} 



■dyp hv d\!/ hz d^. ■ 
The last integrand of i/-,, is 



(90) 



_ / • I A, 3 I j>. + 1 A 23 | yf/kx , I A33 I / | A 13 1 & + 1 A 23 I ^x Vl 
P I J I A 33 I + 2 I A I I I Ass I / J 

X f e; 

t/ — to 



-P\ yp hy + 



A 13 [ \j/ r + 1 A 23 I ^ 



exp - 



A 33 



2 I A 



(*+ ' i -'*it'| i "'*') f ) 



>(ty 



r *« 



. exp U (l *» I *• + K- I *» 
V V I A 33 [ 

I A 33 j / | Ai 3 1^+1 A 23 I \p h:i 
2 I A I \ I A 33 I 



+ 



*£«»(*-M* 
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From Cramer, 15 p. 99 we obtain 

L exp (* ~ t f h = vlf exp ~k ' 

where 

/, _ I A 33 |. 

Then following the same techniques we find 

M *s» = o*t\ a i* £ exp («• ~ 2 *)* 

•/" exp (j^ 2 - |^)^ 2 £ exp (jtf, - |^ 

L 2 W ^ A 2 ^ A 3 /l 



where 



(2tt) j 

A, =■!■ 

/i 2 = 7 
^3 = T 



hihoha 



exp 



B 



A I | A, 

1 



A 



A 33 \B 



[BC - A 2 ] 



Vi = v 



v 2 = 1 



V-, = 



1 - 



A;,; 

A, 



-11- 



A„ \) A 



A = I Ai 2 I I A 33 I — I A i: 



= -pi* I A 



B = | Asa | 

C = | A„ | 
5C - A 2 = 



A33 I I A 2 ; 

A 33 I — I A,j 
I Aas I I A I 2 . 



= p'u I A 

= P22 I A 



The constant value outside the bracket of (91) is 

1 /SL__L_ / 

(2tt) j I A p MhJUi, " I A p "V 



£C - A' 



= 1 



(91) 
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and the expression inside the bracket of (91) is 



2U 



I i <" i v f\ 

exp -o I it + r + r J 



2 \_h, ' h s h a 



1 I A 



, , (I A 33 I — I A 23 I)" 

1 "*" B 



"' Xp l 2|A 33 

I A 33 I - I A l8 I - (| A 33 I - I A 23 |) I] *B 



+ 



= exp {-|(pj, + p^ 2 + P33 + 2pi'> + 2 P f> 2 } = exp -§p{» 2 . 



Thus, 



ilf *,(*>) = exp -|pjw . 



(92) 
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